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We perform fully dynamical simulations at small quark masses by reweighting in the quark mass, 
calculating the weight factors stochastically. This approach avoids some of the technical difficul- 
ties associated with direct simulations. We find that the weight factors fluctuate only moderately 
on nHYP smeared dynamical Wilson-clover ensembles, and demonstrate that the overlap between 
the original and reweighted configurations is large both for short and long distance observables. 
We could successfully reweight 16"*, (LSSfm)"* volume configurations from m^ « 20 MeV to 
mqK,5 MeV, and 24^, (2.77fm)^ configurations from m^ ~ 8MeV to m^j k, 3 MeV quark masses, 
reaching the e-regime. Using the pseudoscalar and axialvector correlators we predict the low 
energy constants E and F and study their volume and mass dependence. 
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1. Introduction 

With efficient simulation techniques, improved lattice actions and increased computational 
resources, essentially all parameter values of lattice QCD, including the point of physical quark 
masses, are accessible to direct simulation today. However, in the small quark mass regime the 
challenges are still considerable: 

• Large volumes are needed for the stability of the algorithms when Wilson fermions are used 

• Autocorrelation times increase dramatically towards the chiral limit 

• Statistical fluctuations of fermionic correlators become difficult to estimate since configu- 
rations with large contributions become rare as small Dirac modes are more and more sup- 
pressed. 

A possible solution to counter these problems is to avoid generating an ensemble with the fermionic 
weight of the light target quark mass but instead simulate a heavier quark mass and reweight to the 
desired ensemble. Simulational algorithms are more efficient at the larger mass and smaller vol- 
umes are sufficient from the algorithmic point of view, while the autocorrelation time is controlled 
by the heavier simulational quark mass. At a larger quark mass the region of small Dirac eigen- 
values is oversampled with respect to the target distribution, thus observables that receive large 
contributions there (e.g. pseudoscalar correlators) will be better estimated. 

On the other hand reweighting will fail if the overlap between the generated and target config- 
uration ensembles is small, or if the weight factors fluctuate strongly. One also has to overcome the 
computational problem of calculating the weight factors effectively. In two recent papers [^, |^ we 
have developed a technique that solves the latter problem and demonstrated that reweighting is pos- 
sible and efficient on fairly large volumes between significantly different quark masses. We have 
applied the technique to reach the £ -regime with Wilson fermions and predicted the values of the 
low energy constants £ and F . In this paper we summarize the technique, giving further supportive 
evidence for the validity of reweighting. We also briefly discuss the £-regime analysis, including 
the use of Wilson fermions in the £-regime and the validity of the continuum chiral perturbation 
theory in our analysis. 

An alternative approach to reweighting was presented in this conference [||]. 

2. Reweighting in the quark mass 

The numerical simulations for this project were done with 2 flavors of nHYP smeared Wilson- 
clover fermions and one-loop Symanzik improved gauge action[Q]. We use the tree-level csw = 1.0 
clover coefficient, so our action is not fully i?{a) improved, but we expect that the 0'{a) corrections 
are small. We have generated two sets of gauge ensembles, both at gauge coupling j8 = 7.2. The 
first set consists 180, 16^ configurations at K" = 0.1278, the second 154, 24^ configurations at 
JC = 0.12805. 

We set the lattice scale from the static quark potential, using ro = 0.49 fm for the Sommer 
parameter. On both configuration sets we found ro/a = 4.25(2) , giving a = 0.1153(5) fm. With 
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Figure 1: The distribution of the plaquette (left panel) and scalar correlator at f = 8 (right panel) on the 
16'* ensemble at K" = 0.1280. On both panels the dashed line is the original (partially quenched) distribution 
and the solid lines correspond to the reweighted distributions. The reweighting included the jSnHYP term as 
explained in the text. 

this value the physical volumes are (1.85 fm)^ and (2.77 fm)^. We connect the PCAC quark mass to 
the renormalized one using the pseudoscalar and axialvector renormalization factors obtained with 
the RI-MOM method, Z'^^ = 0.90(2), = 0.99(2), and estimate the renormalized quark mass 
in the MS scheme at 2 GeV to be 22 and 8.5MeV on the 2 configuration sets, respectively. 

Starting from the original configurations one can explore a range of quark masses in fully 
dynamical systems by reweighting the configurations. The necessary weight factor is the ratio 
of the fermion determinants that can be calculated stochastically, but one must take care not to 
introduce significant statistical errors with the stochastic process. We apply three methods, low 
mode separation, determinant breakup, and ultraviolet (UV) noise reduction to control the statistical 
fluctuations. In both the 16^ and 24^ ensembles we separate 6 low Hermitian eigenmodes. In 
addition we break up the determinant to the product of 33 and 60 terms for each Ak" = 0.0001 shift 
in reweighting on the 16^ and 24^ volumes, respectively. 

To control and remove some of the UV noise we introduce an nHYP plaquette pure gauge term 
in the reweighted action with coefficient jSnHYP = 6.0(k: — JQew)- This value is so small that within 
errors we do not observe its effect on the lattice spacings or quark masses. It is worth emphasizing 
that including the nHYP plaquette gives simply an alternative gauge action and it does not introduce 
any systematical errors. 

With the above outlined method we could reweight the 16^ configurations from K = 0.1278 
to K = 0.1281, or approximately 5MeV quark mass, and the 24^* configurations from k: = 0.12805 
to K" = 0.I28I5, or 3.8MeV quarks. The calculation automatically gives the weight factors at in- 
between K values and we analyzed both data sets at 4 different masses. 

In addition to removing the UV fluctuations and thus reducing the fluctuations of the weight 
factors, the introduction of the nHYP plaquette term also increases the overlap between the orig- 
inal and target ensembles. The largest weight factors are pushed from the edge of the plaquette 
distribution to the middle, where the statistical sampling is better, and the effect is similar for other 
observables as well. We illustrate this in Figure [T] with the distributions of the plaquette and the 
scalar correlator at ? = 8. All data correspond to the 16^ data set at k: = 0.1280 9 MeV quarks). 



3 



Epsilon regime calculations with reweighted clover fermions 



Anna Hasenfratz 



0.03 



o partial quenched 



reweighted 



0.01 — J; 




I - 



0.00 



-0.01 







5 



10 



15 



t 



Figure 2: The scalar correlator atK — 0. 1280 for the partially quenched and for the reweighted ensembles. 

For both quantities the overlap between the reweighted and original distributions is excellent, there 
is no sign that reweighting would prefer region that is poorly sampled by the original ensemble. The 
main difference between the reweighted and partially quenched scalar data is the suppression of the 
negative contribution of the latter one. This seemingly slight difference is nevertheless sufficient to 
make the scalar correlator positive. In Figure ^ we compare the partially quenched and reweighted, 
thus fully dynamical, scalar correlators. Not only does the reweighted correlator stay positive, its 
statistical errors are also significantly reduced compared to the partially quenched one. The overlap 
between the reweighted and partially quenched correlators for other mesons are similar, typically 
better, than for the scalar one. 

3. e -regime analysis 

In the e-regime the pion correlation length is large compared to the linear size of the lattice, 
the light pseudoscalar mesons dominate the dynamics. Nevertheless, in order to incorporate the 
massive modes the volume has to be large compared to the QCD scale, or L » F. Chiral perturba- 
tion theory predicts that at next-to-leading (NLO) order the meson correlators are quadratic in time 
and depend only on two low-energy constants, £ = \imm^o{qq) and F = limm^oFjc- The explicit 
formulas have been calculated for the different meson correlators in continuum xPT, averaged 
over topological sectors, in [§, ^]. In the case of the Wilson lattice action one expects extra terms 
in the chiral expansion, but they show up only at next-to-next-to-leading order, so as long as the 
NLO formulas describe the data, the continuum expressions are sufficient. 

An important question to consider if Wilson fermions can be used at all in the e-regime. As the 
quark mass decreases in a large volume (p-regime) simulation, the chiral symmetry breaking effects 
of Wilson fermions get large compared to the mass, and that can create large lattice artifacts. In 
practice the continuum limit has to be taken before the chiral limit. The situation is different in the 
e-regime, where the finite volume of the system creates an infrared cutoff even at vanishing quark 
mass. This effect is well illustrated by the Hermitian gap distribution. While in infinite volume 
one expects the median of the gap to scale with the mass p. = Z^mpcAC in the e-regime p., 
governed by the IR cutoff of the volume, remains finite while mpcAC — 0. This is clearly the case 
in our simulations. While the quark mass changes by a factor of 4 as we reweight from K = 0.1278 
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Figure 3: The pseudoscalar ( red diamonds) and axialvector (blue bursts) lattice correlators and the com- 
bined fit results on the 16^* data set. The axialvector correlators are multiplied by the factor 50 to better 
match the scale of the pseudoscalar. 



to K" = 0.1281 on the 16^ ensemble, the median of the gap changes by about 40% only One 
does not need a chiral action to study the epsilon regime, though the explicit symmetry breaking 
effects should be small compared to the inverse lattice size. 

In this study we consider the pseudoscalar and the axialvector correlators. The former one is 
dominated by £, the latter by F , though both depend on the &{\) combination mZV and the ^(e^) 
quantity 1 / (FL)^. The lattice correlators have to be multiplied by the renormalization factors Zp 
and Z| to obtain the continuum ones, while in the product mLV the quark mass can be expressed 
in terms of the PC AC mass as m = mpcAC^A/^p- In our analysis we do a combined fit to the 
continuum pseudoscalar and axialvector correlators. 

The £-expansion formulas are systematic expansions in the parameter 1/(FL)^ = iff{e^), but 
depend on the ^(1) quantity mLV. In our simulation we explore the range mLV « 0.7 — 5.0. 
Large values introduce large NLO and NNLO corrections to the correlators, and at some point one 
transitions into the large volume jc-regime. Only by examining the fit results can we decide what 
range of mLV values are acceptable in the £-regime. The volume dependence is also an important 
issue to consider. On our smaller, (l.SVfm)"^ volume, the e-regime expansion parameter is rather 
large, 1/{FL)^ ^ 1.45, and one expects large corrections to the NLO formulas. On our larger 
volume l/(f L)2 0.46, the NLO expansion is much more reliable. 

Our philosophy is to fit both the small and large volume data with the NLO xPT formulas 
and study the quality of the fit and the change of the predicted low energy constants as the volume 
increases and quark mass decreases. 

The fits on the 16^ lattices are shown in Figure ^, where we plot both the pseudoscalar and 
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Figure 4: Same as Figure 3 but on the 24 data set . 
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Figure 5: The low-energy constant F and as the function of the parameter mYV , predicted by NLO 
ZPT. 



axialvector correlators. The data are well described by the NLO formulas at all four mass values. 
Figure |5| shows the predictions for F and £ as the function of mLV . While the condensate £ is 
basically constant, F shows a slight drift as the quark mass decreases. 

Our second data set is 24^, (2.77fm)^, considerably larger. Figure ^ shows the combined fit 
and, again, we find good agreement for all correlators. The statistical errors are under control 
everywhere, though they increase as the reweighting range increases . 

The low energy constant £ on the larger volume is consistent with the value obtained on the 
small volume (see Figure In case of F the large volume data is independent of the quark mass 
and is consistent with the small volume data where the two sets overlap. We cannot exclude the 
possibility that in the large volume data we see the breakdown of the £-expansion at large mLV , but 
we believe finite volume corrections are more likely to explain the difference. Since the next-to- 
leading order corrections to F are over 10% on the 16^ data set, we prefer using the large volume 
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data to arrive at our final prediction, 

F = 90(4)MeV, = 248(6)MeV . (3. 1) 

Tlie errors only include the statistical uncertainties. These values are consistent with other recent 
two-flavor computationsi[^, 



4. Conclusion 



The data presented in this paper have been generated with moderate computer resources. This 
was possible due to the good chiral properties of the action which come at relatively low cost due 
to the simple nHYP smearing procedure, and the effective reweighting that allowed us to lower the 
quark mass even further. We have demonstrated good overlap for different observables between 
the original and reweighted data sets. 

With reweighting we were able to reach the £ -regime with Wilson fermions on 2 different 
volumes. Since the ;^PT expansion converges slowly in 1/(FL)^, our large volume puts us into a 
good position and the comparison between the L/a = 16 and L/a = 24 results shows that the finite 
volume effects are under control. 

For the e expansion to be valid, the parameter miy has to be i^(l). Our data span the range 0.7 
to 5.2 and might go beyond the validity of the analytical expressions. An expansion that connects 
the £ and p regimes would be very useful to control this aspect of the calculation. 

Repeating the calculation at a smaller lattice spacing would not be prohibitively expensive and 
could improve on all of the above mentioned issues. 
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